Pion photoproduction off the proton in a gauge-invariant chiral unitary framework 
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We investigate pion photoproduction off the proton in a manifestly gauge- invariant chiral unitary 
extension of chiral perturbation theory. In a first step, we consider meson-baryon scattering taking 
into account all next-to-leading order contact interactions. The resulting low-energy constants are 
determined by a fit to s-wave pion-nucleon scattering and the low-energy data for the reaction 
n~p — > r\n. To assess the theoretical uncertainty, we perform two different fit strategies. Having 
determined the low-energy constants, we then analyse the data on the s-wave multipole amplitudes 
Ea+ of pion and eta photoproduction. These are parameter-free predictions, as the two new low- 
energy constants are determined by the neutron and proton magnetic moments. 
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I. INTRODUCTION AND SUMMARY 



Pseudoscalar meson photoproduction off protons is one 
of the premier tools to unravel the spectrum and prop- 
erties of baryons made of the light up, down and strange 
quarks, as witnessed e.g. by the dedicated baryon reso- 
nance programs at ELSA (Bonn) and CEBAF (Jefferson 
Laboratory). Some of the low- lying resonances like the 
Roper iV*(1440), the 5u(1535) or the A(1405) exhibit 
features that can not easily be reconciled with a simple 
constituent quark model picture. Therefore, it was spec- 
ulated since long that some of these peculiar states and 
their properties can be explained if one assumes that they 
are generated through strong coupled-channel dynamics. 
Arguably the best tool to address such a dynamical gen- 
eration of resonances is unitarized chiral perturbation 
theory [H-i]. 

Over the years, we have developed and a applied a 
gauge-invariant chiral unitary coupled-channel approach 
based on the leading order (LO) chiral effective La- 
grangian of QCD to kaon [l(| and eta photoproduction 
[ll|. To go beyond LO, one first has to refine the de- 
scription of meson-baryon scattering in this framework as 
the strong hadronic final-state interactions are a crucial 
ingredient in evaluating the complete photoproduction 
amplitudes. Therefore, in a recent work Q we have de- 
veloped a framework to analyse meson-baryon scattering 
incorporating next-to- leading order (NLO) contributions 
of the chiral Lagrangian Our scheme is based on 

the solution of the Bethe-Salpeter equation (BSE), with 
a kernel derived from the contact terms of the NLO chiral 
Lagrangian. We have summed up the full infinite series of 
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Feynman diagrams generated by the BSE without resort- 
ing to any of the commonly made approximations as e.g. 
the on-shell approximation. In this way we were able to 
reproduce successfully both s- waves of pion-nucleon scat- 
tering (Sn and 531) between the irN and r]N thresholds. 
For higher energies we have observed that only the reso- 
nance Sn (1535) but not the S31 (1620) could be described 
well. We have concluded that the S3i(1620) does not 
have a prominent dynamically generated content. As a 
matter of fact after fixing the Sn partial wave in the first 
energy region, i.e. W cms < 1.56 GeV, the Sn amplitude 
for the higher energies came out in astonishing agreement 
with current partial wave analyses. The examination of 
the complex energy plane showed that also the pole po- 
sition of the second s-wave resonance, i.e. the 5*11(1650), 
agrees rather well with those given by the particle data 
group [l3| . It is therefore natural to extend this approach 
to s-wave photoproduction of pions, where a large data 
basis already exists. This is done in this paper. 

The main results of our investigation can be summa- 
rized as follows: 

1) As the first step, we have considered the meson- 
baryon scattering processes that are relevant for 
pion photoproduction. We consider six coupled 
channels, utilizing the chiral effective Lagrangian 
at NLO. In the kernel of the underlying Bethe- 
Salpeter equation, we include all terms allowed by 
the symmetries. We do not use the on-shell ap- 
proximation that is common in most works on re- 
lated subjects. The crossed-channel contributions 
are not yet included since the proper treatment of 
the corresponding diagrams in a framework com- 
parable to the one presented here constitutes an 
unsolved problem. 

2) To pin down the parameters of the approach (low- 
energy and subtraction constants), we perform two 
fit strategies. In strategy (I), we use as input the 
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data on elastic irN scattering in the Sn and Ssi 
partial waves for energies in the range {rrip + M^) < 
Wcms < 1.56 GcV. In strategy (II), three subtrac- 
tion constants are fixed and the Sn partial wave is 
fitted up to W cms < 1.70 GeV, but the 5 3 i only in 
the near-threshold region, W cms < 1.20 GeV. The 
data on ir~p — > r\n from Rcf. flij is included in 
both fit strategies. 



3) In both fit strategies, the S\\ partial wave and 
the data on n~p — > r/n are well described. The 
7V*(1535) and the iV*(1650) are both dynamically 
generated, the precise pole positions depend on 
the fit strategy, cf. Eqs. (|11I14[) . We also give 
predictions for the scattering lengths a n N^riN and 

4) Having scrutinized the hadronic sector, we have 
extended our approach to s-wave pion photo pro- 
duction. The only new parameters can be deter- 
mined from the nucleon magnetic moments and 
thus parameter-free predictions emerge. We find 
a good description of the s-wave multipolc £'o+ for 
pion photoproduction in the .Sn-wave and also for 
r] photoproduction. 

Having summarized the most important results of our 
study, it is important to briefly discuss possible improve- 
ments of the method. First, the crossed channel dynam- 
ics has to be included properly. This will allow e.g. to 
get a better description of the near-threshold region in 
pion- nucleon scattering. Unfortunately, the exact imple- 
mentation of both crossing symmetry and unitary has 
not been possible so far in approaches based on Feyn- 
man diagrams, in contrast to other approaches as e.g. 
the one based on Roy-Steiner equations [15j where these 
constraints are met by construction. For an attempt to 
approximately restore crossing symmetry in an ansatz 
comparable to the one employed here, see e.g. [l^ . Fur- 
ther, in some channels, explicit resonance degrees of free- 
dom will have to be incorporated as not all resonances 
are generated dynamically. For a method to do that, 
see e.g. Ref. [131 ■ Finally, a larger data base including 
also kaon-nucleon scattering and kaon photoproduction 
should be considered simultaneously with the processes 
studied here. Work along these lines is in progress. 

The paper is organized as follows: In Sec. [Til t nG un- 
derlying approach to analyse meson-baryon scattering is 
described. Two fit strategies to pin down the occurring 
parameters are presented and fits to data and predictions 
for some scattering lengths and the fundamental rjN scat- 
tering amplitude are given. The approach is extended to 
photoproduction in Sec. IIII1 where parameter-free pre- 
dictions for pion and 7/ production off protons are given. 
Many technicalities are relegated to the appendices. 



II. HADRONIC SCATTERING 

In the present section we wish to describe our frame- 
work for meson-baryon scattering as it will serve us as the 
most important ingredient in our analysis of the photo- 
production processes. We wish also to give a more sys- 
tematic error analysis of our results as it was done in [7| ■ 
Moreover from this first study we know that our coupled 
channel approach is applicable to relatively high ener- 
gies, thus we will extend our analysis also to the process 
ttN rjN. 

A. Formalism 

Chiral Perturbation Theory (ChPT) is a powerful tool 
to analyse pion-nucleon scattering and pion photopro- 
duction at low energies in a systematic manner. In 
its original formulation one starts from the chiral La- 
grangian. The reaction amplitudes are then given by a 
set of diagrams including contact interactions as well as 
loop diagrams. Working order by order in a low-energy 
expansion, every Feynman graph can be classified by the 
order to which it contributes in this expansion. For a spe- 
cific chiral order the number of such diagrams is finite. 
Thus technical effort for a calculation of the scattering 
amplitudes can be large but remains finite. Up to now the 
scattering amplitudes for meson-baryon scattering have 
been calculated in the three-flavour formulation of ChPT 
up to the third chiral order in Refs. [H, [l||. However, 
the region where the result of this strictly perturbative 
approach yields a decent description of the process is re- 
stricted to the subthreshold region, and it certainly fails 
in the vicinity of resonances (as nicely discussed e.g. in 
Ref. 0). 

In the following we present a framework to analyse 
meson-baryon scattering for energies up to the resonance 
region. It is inspired by Chiral Perturbation Theory, 
which gives the driving terms of the meson-baryon in- 
teraction, and it relies on the coupled channel Bethe- 
Salpeter equation (2lj . The latter implements the re- 
quirement of two-body unitarity exactly and in principle 
allows to generate resonances dynamically. It already 
improved our understanding of the purely mesonic and 
meson-baryon sector in recent years, see e.g. 

We denote the in- and out-going meson momenta 
by qi and q 2 , respectively. Moreover the overall four- 
momentum is given by p = qi + p\ = <?2 + P2, where 
Pi and P2 are the momenta of in- and out-going baryon, 
respectively. For the meson-baryon scattering amplitude 
T(q2,qi;p) and potential V(q2,qi;p) the integral equa- 
tion to solve reads in d dimensions 

T(q2,qi;p) = V(q 2 ,qv,p) (1) 

+ i J 4Jy[V(q 2 ,l-,p)Sty- t)A(l)T(l, qi ;p), 

where S and A represent the baryon (of mass m) and 
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the meson (of mass M) propagator, respectively, and 
are given by iS(p) = i / — m + ie) and iA(fc) = 
i/(k 2 — M 2 + ie). Let us note that all objects in the 
above equation which are denoted by capital roman 
letters, are matrices in the channel space. Restrict- 
ing the meson-baryon states to have the quantum num- 
bers of the proton, the channel space reduces to the 
6-dimcnsional space spanned by the following channels 
{TT°p,iT + n,rip, A"+A, A' + £°, A"°£ + }. As a matter of fact, 
the BSE can be visualized as presented in FigJTJ Note 
that the loop integration in Eq. ([I} is performed in d di- 
mensions, without restricting the loop momenta to be on 
the mass shell. Such an approximation would certainly 
reduce the technical effort to solve the BSE, however it 
spoils the direct correspondence of the solution of the 
BSE Eq. ([1]) to the series of Feynman graphs, which we 
evaluate as it stands. Thus, every term in our solution 
of the BSE is related, in a one-to-one correspondence, to 
a properly evaluated Feynman diagram. It comes as an 
advantage of this prescription that the implementation 
of gauge invariance in a combined analysis of hadronic 
scattering and meson photoproduction is straightforward 
and very natural: It follows essentially the guidelines 
from quantum field theory textbooks, see e.g sec. (7.4) of 
[22j , and will be explained in sec. IIII1 Moreover, it is also 
straightforward to compare our amplitude to the pertur- 
bative expansion at any fixed chiral order. We complete 
our discussion of the off-shell-dependence of the effective 
vertices with the remark that the use of the on-shell- 
approximation is not more "physical" than taking the 
off-shcll-dcpendcnce into account, though reducing the 
effort of the calculation significantly down to the evalu- 
ation of a geometric series. Simply iterating a fixed on- 
shell kernel in such a geometric series can even lead to 
significant deviations from the results of Feynman graphs 
when iterating Born-terms, as is exemplified by an anal- 
ysis of box graphs in sec. (5.2) of In our case, off- 
shell behavior of the potential reflects itself in tadpole 
- integral terms in the full scattering amplitude. These 
terms might in general depend on the chosen parame- 
terization of fields as long as only a subset of Feynman 
graphs is summed up. Setting them to zero as done us- 
ing the on-shell approximation is just one possible "choice 
of gauge" (in the space of field parameterizations) in a 
non-invariant result, which however is not in line with 
the proper evaluation of loop diagrams we aim at here. 
As the analytic energy-dependence of the tadpole-terms 
is trivial, it should be possible to compensate for this 
non-invariance effect by an adjustment of coupling con- 
stants in the kernel. As a byproduct of our procedure 
for finding fits 1 , we have checked numerically that this is 
indeed the case. Thus, while the exact numerical values 
of the coupling constants should be taken with a grain 
of salt (they should be considered as model parameters 



FIG. 1: Symbolical representation of the Bethe-Salpeter equa- 
tion. Here the circle and the square represent the potential 
V and the scattering amplitude T, respectively. 



in our approach) the overall properties of the amplitude 
are solely based on the unambiguous analytic properties 
of the selected infinite subset of loop graphs. 



B. Interaction kernel 

The only missing part of equation (TI} remains the driv- 
ing term of the meson-baryon interaction, i.e. the poten- 
tial V(q2, <7i;p). It is well known that already the leading 
order chiral potential, the so-called Weinberg- Tomozawa 
term, captures the prominent part of the meson-baryon 
dynamics at the ttN threshold. It is derived from the co- 
variant derivative of the leading order chiral Lagrangian, 
which reads 



4b = (Bii-f.D* - mo)B) 



D/F 



(B lfll5 [u»,B}±), 



(2) 



where (...) denotes the trace in flavor space, [D^, B] := 
d^B + ^[[u^, d^u], B], too is the common baryon octet 
mass in the chiral limit, while D and F are the lowest- 
order axial coupling constants. The relevant degrees of 
freedom of ChPT are the Goldstone bosons, which are 
collected in the traceless meson matrix U <G SU(3), 
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where Fo is the meson decay constant in the chiral limit. 
The baryonic fields are collected in a traceless matrix 
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EL 4- A. 
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In view of the application to meson-baryon scattering we 
set external currents to zero except for the scalar one, 
which is set equal to the quark mass matrix, s = M. := 
diag(TO„, TO(2, to s ). Note that this situation will change 
as soon as we will deal with photoproduction in the next 
section. We furthermore use 



1 See the discussion in sec. 3 of Q 



u 2 :=U, u^ 1 := iv)d^u - iud^u\ 



f ± ux V X 



2B s 



(5) 
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where the constant Bq is related to the quark condensate 
in the chiral limit. Separating the channel space from the 
momentum space structures, the Weinberg- Tomozawa 
potential reads 



V W T{q2,qi;p) 



(6) 



where Awt denotes a matrix in channel space as de- 
fined in App. [AJ The axial coupling in the Lagrangian 
([2|) induces an additional contribution to the chiral po- 
tential formally at the leading chiral order. It arises via 
the one-baryon exchange graphs in the u and s-channel, 
the so-called Born graphs. It is not known how to fully 
include such terms in the scattering kernel of the BSE 
when attempting an analytic solution in d dimensions, 
corresponding to Fcynman graphs. The difficulties one 
encounters in this case are described in Q . On the other 
hand, it is plausible that the s-wave contributions stem- 
ming from the Born terms can be well approximated by a 
series of contact terms not too far above the two-particle 
thresholds. We shall return to this open problem in a 
separate publication [24[. In the present work, wc will 
restrict the driving term of our coupled channel analysis 
to a set of contact terms. We go beyond the leading or- 
der and complete this set by the meson-baryon vertices 
from the second order chiral Lagrangian from (25l - [27| . It 
reads 

cf B = bo(BB)( x+ ) + b D/F (B[ x+ ,B] ± ) + b 4 (BB)( Ufl u») 
+ 6i/ 2 (b[«/*, [«", B] T ] > + b 3 (B{ Ufl , {< B}}) 



+ tb 5/6 (B<T^ 



[ Ufl ,u v ],B ) +ib 7 (Ba* iU Uu)(u v B) 



^((^[^,K,[^,B]] T ]) 



(B-f D v ,[u u ,[ Ufl ,B]} : 



+ ^((^K,K,[^,b]}}) 



(By[u„K 1 K 1 B}}])) 



^(2(B^[D v ,B])(u^) 



(B-fB) { [D v , u„] u v + u M [A, , «"] ) 



M^[/r>s]>+M^{/r>s}>), (7) 



where fit" includes the electromagnetic field strength 
tensor, which vanishes for u M = but will become impor- 
tant for the photoproduction amplitude later. All opera- 
tors of the second chiral order in Eq. (0) are accompanied 
by the low-energy constants (LECs) bi. We complete the 
driving term Eq. © of Eq. (JTJ by the second order chiral 
potential, which reads 

VNLo(q2,qi;p) =An(q 1 ■ q 2 ) + A 57 [<fi,qk] + A M 

+ A 811 (q/ 2 (q 1 -p) + fi(q 2 -p)), (8) 



where the matrices A... in the channel space depend on 
the LECs as detailed in App. [A] As a matter of fact 
the importance of the second order terms in the kernel 
of the BSE is twofold. First of all, as can be seen in 
[lH j such terms lead to sizable corrections of the meson- 
baryon scattering amplitudes. Secondly, the contact in- 
teractions of the second chiral order not only contribute 
to the s-wave but also to the p- waves, which are then 
iterated in the BSE. In App. [D] we demonstrate in a toy 
model that the presence of the first two partial waves is 
sufficient to reproduce the correct behavior of the differ- 
ential cross sections at sufficiently low energies. 



C. Solution of the BSE 

For the solution of the BSE we use the recipe devel- 
oped in a foregoing publication Q. As described there 
we utilize dimensional regularization to treat the diver- 
gent loop integrals, where the purely baryonic integrals 
are set to zero from the beginning, while only an energy- 
independent constant is subtracted from the fundamental 
meson-baryon loop integral. This prescription to treat 
the large baryon mass scale is similar to the EOMS reg- 
ularization scheme described in p8j . On the other hand 
it resolves several technical problems, which appear in 
the course of the study of the photoproduction ampli- 
tudes described in the next section. To further extent, 
the solution of the BSE corresponds to an infinite chain of 
one-meson-one-baryon loop diagrams, see Fig. [TJ From 
the point of view of the usual perturbative treatment this 
would demand an infinite number of counterterms from 
a local Lagrangian to absorb the loop divergences. This 
is of course not feasible in an effective field theory. In 
our non-perturbative framework, the ignorance of higher- 
order terms in the scattering kernel, which would serve 
to cancel the divergences and the scale-dependence of the 
loop integrals in a perturbative setting, reflects itself by 
the appearance of a new free parameter for every loop- 
integration, parametrized here by the logarithm of the 
rcnormalization scale. This pragmatic approach is com- 
monly adopted in the literature, see e.g. @-H, H3|- 
The new free parameters arc not completely arbitrary, 
however: At least, we must impose that the values for 
the renormalization scale correspond to neglected higher 
order terms of natural size. Should this not be the case, 
and a scale of e.g. \i ~ TeV emerge from some fits, we 
must discard that solution as unnatural. As a side re- 
mark, we note that any modification of the loop integrals 
corresponds to a specific modification of the potential V 
in the solution of the BSE. For an explicit demonstra- 
tion of this procedure we refer to App. F of [3l|. The 
requirement that the modification of the potential is not 
dominating the leading order terms also yields the men- 
tioned constraints on the free scales. In conclusion, the 
foregoing discussion suggests that it is sufficient in the 
present work to apply the subtraction scheme described 
above, keeping in mind that the modified loop integrals 
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still depend on the renormalization scale, which consti- 
tutes a free fitting parameter. In the next subsection we 
will re-examine this as well as the possibility to adjust 
this scale to a fixed value due to constraints on the loop 
dressing of vertex functions. 

The essential advantage of the above treatment is the 
preservation of the analytic structure of the loop inte- 
grals, which allows for a continuation of the scattering 
amplitudes into the complex energy plane. The solution 
of the BSE is presented in App. [5] It can be written in 
terms of elementary functions (that is without resorting 
to a numerical solution) of the loop integrals, which are 
collected in App. [Cj Once the BSE has been solved we 
put the external particles on their mass shell and calcu- 
late partial wave amplitudes as well as differential cross 
section for hadronic scattering. For the evaluation of 
the photoproduction amplitudes we will require the full 
off-shell dependence of the hadronic solution, as will be 
described later. 



D. Fit 

It is important to clarify the physical input to the scat- 
tering amplitudes. Throughout this work we will use 
physical hadron masses (in GcV), i.e. M^o = 0.1350, 
M w+ = 0.1396, M v = 0.5478, M K + = 0.4937, M K o = 
0.4977, m p = 0.9383, m„ = 0.9396, m A = 1.1157, 
mjo = 1.1926 and m s + = 1.1894. The baryon mass in 
the chiral limit, mo in Eq. ([7]), can be fixed to 1 GeV with- 
out loss of generality, as any other value only amounts 
to a rescalin g of the unknown LECs. In contrast to the 
earlier work [lOfl , the meson decay constants are fixed to 
their physical values, i.e. F n = F n /1.3 = 0.0924 GcV, 
F K =0.113 GcV. 

To pin down the free parameters of our approach we 
have to specify experimental input available on the mar- 
ket for the considered meson-baryon channels. From the 
experimental point of view elastic ttN scattering is by far 
the best explored reaction. On the other hand it is clear 
that the low-energy region is dominated by the p-wave 
resonances, namely Roper and Delta. Our foregoing in- 
vestigations in [3] have shown that we are not able to 
dynamically generate these resonances consistently with 
the s-wave resonances. Since those degrees of freedom 
are not included in our approach, we restrict ourselves to 
the analysis of s- waves. We fit our results for these irN 
partial waves to the widely accepted partial wave anal- 
ysis (WI08) by the SAID collaboration [33. Comparing 
an earlier analysis by the Karlsruhe group |33j to the cur- 
rent one, we assign for the energies below W — 1.28 GcV 
an absolute systematic error of 0.005 and for higher en- 
ergies an error of 0.030 to the partial wave amplitudes. 
To some extent this is in agreement with error estimates 
done in [29|, which are motivated by the expectation of 
pronounced three-body effects above the irirN threshold. 

Another widely explored channel is 7r _ p — > rjn, for 
which we consider quite recent but already very estab- 



lished results on differential cross sections measured by 
Prakhov et al. in [l4j]. For all seven measured incident 
pion beam momenta p\ a b we assign a measurement er- 
ror as well as the systematic error of 6% as pointed out 
in Ref . [l4[ ■ Moreover one should keep in mind that also 
Piab itself entails an uncertainty, which hampers the clear 
pairwise separation between most of the given beam mo- 
menta [34| . We do not consider this uncertainty in our 
fitting routine as the inclusion would require an addi- 
tional model-dependent input. The necessary formalism 
is collected in App. [Dl 

To fit the above data we follow two different fit strate- 
gics. Fit strategy (I): We start from the best fit ob- 
tained in and additionally include the n~p — > rjn dif- 
ferential cross sections by Prakhov et al. adjusting all 
17 parameters of the model (log(fi 7r ), log(^_R-), log(/i I( ) as 
well as 14 LECs). The fitting region in the elastic irN 
channel is chosen to be (m p + M n ) < W cms < 1.56 GeV 
for both Sn and S31. It is obvious that the new data 
will restrict our parameter set additionally, possibly cor- 
rupting the agreement of the elastic irN to the SAID 
data compared to the fit obtained in m. Fit strategy 
(II): One of the main observations of [7j was that the 
Sn but not the 1S31 partial wave of 7r7V scattering can 
be described well in the resonance region. On the other 
hand the main goal of the present work is to see how 
hadronic resonances manifest themselves in the photo- 
production amplitude. Moreover since the rjN final state 
is an isospin 1/2 state we fit the elastic Su partial wave 
in the energy region {m p + M„) < W cnls < 1.7 GeV 
together with ir~p —¥ r\n differential cross sections by 
Prakhov et al.. The 6*31 is considered only in the near- 
threshold region for W cms < 1.2 GeV. Furthermore, we 
reduce the number of the free parameters of the model 
to LECs fixing the regularization scales from the begin- 
ning. It turns out that the corrections to the tree-level 
result of the photoproduction multipole i?o+ due to the 
dressing of the B — > <f)B vertex are large already at the 
7riV threshold. Thus we choose the scales such that the 
meson-baryon loop integral evaluated at s = nip vanishes 
in every meson-baryon channel, to assure that the axial 
vertex-function does not deviate much from the corre- 
sponding tree-level expression. 

For both fit strategies we minimize the following quan- 
tity 

2 2 J2i n i fn\ 

x -=x DOF = N{Eini _ p) i:-- o) 

Here p is the number of free parameters, rii is the num- 
ber of data points available for the observable i and TV is 
the number of observables. This choice of \ 2 ensures the 
equal weight of both fitted observables, compensating for 
the different number of data points. Albeit the solution 
of the BSE is fully analytical, it costs a huge amount of 
computational power. Thus the fitting as well as error 
estimation routine is performed in a parallelized version 
on 20-30 threads of the HISKP cluster utilizing the (mi- 
grad) minimization routine of the MINUIT C++ library 
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|35| . The uncertainty of the model is estimated as follows. 
First, after obtaining the best fit (Xbest) t ne errors 011 
the parameters are calculated in the (hesse) subroutine of 
the MINUIT package. Then within these errors we gen- 
erate a large number of parameter sets (~ 10, 000) and 
calculate for each the corresponding Xdof- Then each 
set corresponding to a \ 2 < Xbest + 1.15 is considered 
to produce results in the la region around the central 
value 2 . 



E. Results 

Solution I: Following the first fit strategy, we obtain 
the best fit as presented in Figs [3] As already dis- 
cussed, the differential cross sections on n~p — > i]n are 
observed at seven different pion beam energies, which by 
themselves entail a non-negligible uncertainty. The lat- 
ter, however is not included into the definition of the x 2 
for the reasons given above. Therefore, we refrain from 
giving any numerical value for this quantity 3 . The cor- 
responding parameters (bi in GeV -1 and jii in GeV) are 
of natural size and read: 



log(At ff /(lGcV)) = +1.003 ±0.331 
log(/i, ; /(lGeV)) = +1.034 ±0.298 
log(/i K /(lGeV)) =-0.168 + 0.080 



h 


= -0.126 + 0.039 


b 8 = 


+0.610+0.012 


h 


= -0.139+0.045 


b 9 = 


-0.677+0.037 


b 3 


= -2.227+ 0.133 


bio = 


+2.027 + 0.100 




= -0.288+ 0.080 


6n - 


-0.847 + 0.027 


h 


= -1.402 + 0.094 


h = 


-1.063+0.038 


h 


= +0.474+0.118 


b D = 


+0.771 + 0.042 


b 7 


= -1.751 + 0.368 


b F = 


-0.169 + 0.054 



The observation in the elastic irN channels is similar 
the one made in Between the nN and r/ N thresh- 
olds both partial waves can be fitted nicely to the SAID 
partial wave analysis. For the ttN scattering lengths of 
isospin /, etj (in units of 10 _3 /M 7r +), we obtain: 

a^ = -87.0it| and a\% = +174.5™. (10) 

The theoretically cleanest determination of these observ- 
ables stems from the analysis of pionic hydrogen and pi- 
onic deuterium data based on effective field theory [36| . 
Oi/2 = (179.9 ± 3.6) x 10~ 3 /M W+ and a 3/2 = (-78.5 ± 
3.2) x 10 _3 /M 7r +, which is in nice agreement with our de- 
termination for the 7=1/2 channel, but is slightly too 
small for I = 3/2. For both isospins our determination 



2 One might argue, whether or not one should divide by the num- 
ber of degrees of freedom. It turns out that the error bars do not 
change significantly. 

3 We wish to note that \ 2 restricted to the SAID data lies only 
slightly above the one given in the foregoing publication [7|, 
where no other than elastic ttN channels were included as ob- 
scrvablcs. 



agrees perfectly with those from the direct extraction by 
the SAID collaboration: a 1/2 = (174.7+2.2) x 10" 3 /M T + 
and 03/2 = (-89.4+1.7) x lQ- 3 /M n +. 4 

In the higher energy region the lowest Su but not S31 
resonances could be reproduced as dynamically generated 
states in our model. The pole position can be extracted 
via analytic continuation of the scattering amplitude to 
the complex W cms plane and read 

Wl535 = (1.547t°;°j* - z0.046±° a ™ 4 7 ) GeV, 

W 1650 = (1.597t°;°2o - i0.045l° :ul s) GeV. (11) 

Obviously the inclusion of rjN observables into the analy- 
sis forces the pole of AT* (1535) to the higher and the pole 
of A* (1650) to the lower energies compared with the pre- 
vious analysis • This observation is in some agreement 
with the analysis in (37j . There in a meson-exchange 
model the analysis of the inelasticities has shown that a 
simultaneous description of the 77A and ttN scattering 
amplitude is hampered by the missing ttttN channels, 
which arc also missing in our approach. However let us 
repeat that the starting values of the present fit strategy 
are chosen to be those from Q. Although there is no 
reason to doubt about them for the elastic ttN scattering 
one should keep in mind that an inclusion of additional, 
i.e. rjN data alter the amplitudes in the irN channel as 
well. 

For the pion induced eta production, Fig. [3l we ob- 
serve that the outcome of the model agrees with the ex- 
perimental data, keeping in mind the uncertainty on the 
pion beam momenta. The inclusion of the latter is cru- 
cial especially for the lowest pion beam momenta, where 
the slope of the total cross section is enormous as can 
be seen in Fig. 0] There we present the outcome of the 
model for higher beam momenta than included in the fit. 
Obviously the total cross section agrees with the exper- 
imental data within the error bars. We also observe a 
large qualitative agreement of the outcome of the model 
with the older and less precise measurements, selected 
by the SAID collaboration. The overshooting of the to- 
tal ttN — > r/N cross section by 30% as discussed in [37j 
is overcome obviously via moving of the A* (1535) peak 
towards higher energies, which is a direct consequence 
of the present fit strategy. In App. [D] we show that the 
present model is in principle capable to simulate a cos 2 9 
like behavior in the differential cross sections, usually re- 
ferred to as the influence of the -D13 resonance. The ob- 
servation to be made from Fig. [3] is that the inclusion 
of the elastic ttN channels prevents (or at least damps) 
such behavior. One should notice that the curvature in 
the data is of comparable size as the uncertainty in the 
data as well as the uncertainty band of our approach. 

We can make a prediction of the s-wavc amplitudes for 
clastic rjN as well as ir~p — > r\n scattering as presented 



4 We thank Ron Workman for providing us with these values. 
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FIG. 2: Best fit according to fit strategy (I) to the real and imaginary part of the Su partial wave compared to the WI08 
analysis done by the SAID collaboration [331 - The dashed vertical lines correspond to the two particle thresholds and the 
bold vertical line limits the energy range, up to which the fit has been performed. The blue and red bands represent the la 
uncertainty of our approach as described in the main text. 
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FIG. 3: Best fit according to fit strategy (I) to the differential cross sections for ir p — > -qn from Ref. [Tj- The error bars of 
experimental data include the systematic error of 6% as arg ued in Ref. QJ]. The green bands represent the la uncertainty of 
our approach as described in the main text. 
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FIG. 4: Total cross section of the best fit for ir~p — > r\n 
according to fit strategy (I). The model is fixed to reproduce 
the differential cross section and thus the total cross section as 
measured by Prakhov et al. (black symbols) . The black curve 
including the uncertainty band is the result of our model. The 
blue circles correspond to older measurements as selected by 
the SAID collaboration which are presented for completeness. 



in Fig. [5l Here and in the future we use the Hohler par- 
tial waves denoted by /o+ in contrast to the S\\ partial 
wave used by SAID collaboration, see e.g. Fig. [2j which 
is the /o+ for isospin 1/2 multiplied by q cms . In both 
channels the real and imaginary part shows a similar be- 
havior to the one shown in [38j . However the position 
of the 5*11 peak is systematically shifted to higher en- 
ergies in our approach, which is again a feature of the 
present fit strategy. For the scattering lengths we obtain 
the following values 



a 



(+0.219 
(-0.234 



+0.047 
0.068 
+0.020 
0.024 



i0.235±°;S£l) fm, 
- Hot) fm- 



i0.129: 



(12) 



There is a large spread in the results on ?7TV scattering 
lengths debated for a long time, see [H] for a nice collec- 
tion of those. One can note that most models predict a 
positive real and imaginary part of the scattering length, 
which is in agreement with our result as well. 

Solution II: As argued above, the number of free 
parameters is reduced in this fit strategy by the three 
rcgularization scales. They are fixed such that for each 

meson-baryon channel (i): lMB(jn^,mi,Mi) = in the 
nomenclature of App. [C] which yields the following val- 
ues 

log(/x 7r ) = -0.368, log(/i ?7 ) = 0.056, \og{n K ) = 0.210 . 

The best fit of the 14 LECs, which are the only free 
parameters of the model in the present fit strategy, is 
presented in Fig. [6] and [7] All parameters are of perfect 
natural size and read including the error bars (in GeV -1 ): 
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FIG. 5: The result of our model in the fit strategy (I) for the 
real and imaginary part of the s-wave scattering amplitude of 
the 7]N (top) and n~p — s- r\n (bottom) channels. The error 
bands represent the uncertainty due to the variation of model 
parameters as described in the text. 



61 


= -0.014 ±0.023 


b s = 


+0.272 + 0.015 


60 


= -0.207 ±0.051 


69 = 


-0.483 + 0.032 


b. 


= -1.063 ±0.032 


bw = 


+1.054+0.021 


64 


= -1.312 ±0.023 


hi = 


+0.328+ 0.015 


h 


= -0.628 ±0.060 


bo = 


-1.228 + 0.005 


be, 


= +0.508 ±0.045 


bo = 


+1.097 + 0.011 


b 7 


= +1.041 + 0.191 


b F = 


-0.858 + 0.011 



In the elastic irN channel the Su partial wave agrees al- 
most perfectly in the whole energy range with the one 
from the analysis by the SAID collaboration. The corre- 
sponding scattering lengths are extracted to be (in units 
of 1Q- S /M W +) 



ttN 
a 3/2 



-93.0^ 



and 



ttN 
a l/2 



-168.9±gl|. 



(13) 



A comparison with the result of other calculations, given 
before, shows the same pattern as in the previous fit. 
Both scattering lengths agree within the error bars with 
the direct extraction by the SAID collaboration and are 
smaller than the values extracted in Ref. [36| . 

Both TV* (1535) and TV* (1650) are reproduced as dy- 
namically generated states of the lowest meson and 
baryon octet states. The pole positions of both TV* reso- 
nances read 



W1535 

Wi 650 



(1.512 ±?-i0.070 t 9 5 ) GeV 



+8 
-7 

(1.715lg 



*0.116±||) GeV. 



(14) 



As a matter of fact we expect the pole positions from the 
present fit strategy to be even more realistic than those 
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FIG. 6: Best fit according to fit strategy II to the real and imaginary part of the Su partial wave compared to the WI08 
analysis done by the SAID collaboration [32|. The dashed vertical lines correspond to the two particle thresholds and the 
bold vertical line limits the energy range, up to which the fit has been performed. The blue and red bands represent the la 
uncertainty of our approach as described in the text. 
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FIG. 7: Best fit according to fit strategy II to the differential cross sections for tv p — s> rjn from Ref. [13 ]. The error bars of 
experimental data include the systematic error of 6% as arg ued in Ref. [3. The green bands represent the la uncertainty of 
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FIG. 8: Total cross sections of the best fit for ir~p — > ryn ac- 
cording to fit strategy (II). The model is fixed to reproduce 
the differential cross section and thus also the total cross sec- 
tion as measured by Prakhov et al. (black symbol). The 
black curve including the uncertainty band is the outcome of 
the model. The blue circles correspond to the older measure- 
ments as selected by the SAID collaboration which are only 
presented for completeness. 



from the previous fit strategy as well as from the anal- 
ysis done in Q, where no physical information was in- 
cluded for energies in the region of the second resonance. 
The pole position of the iV*(1535) is perfectly within 
the uncertainty band presented in [l3j . i.e. W1535 = 
(1.490.. .1.530) - i(0.045... 0.125) GeV. On the other hand 
the position of TV* (1650) differs slightly from the one 
given there, i.e. W 1650 = (1.640. .. 1.670) -i(0.050... 0.085) 
GeV. Note that both bands in [l3[ are mostly based on a 
selection of partial wave analyses. The pole positions 
from two comparable theoretical works read VF1535 = 
1.519 - «0.064 GeV and W 1650 = 1.669 - i0.068 GeV 
from [H as well as W1535 = 1-496 - iO.041 GeV and 
W 1650 = 1.686 - i0.096 GeV from [IJ. 

For the pion induced eta production we observe that, 
taking into account the uncertainty of the pion beam en- 
ergy, all seven differential cross sections agree with the 
data by Prakhov et al., see Fig. [7l Again the cos 2 9 be- 
havior does not appear. We have discussed in App. [D] 
that in principle such a behavior could be reproduced in 
our amplitude by means of enhanced contributions from 
the p-waves, which are iterated in our approach. We 
conclude from this observation that such a behavior is 
excluded in this combined nN and r/N fit. In Fig. [5] wc 
present the total cross section for the same process be- 
yond the fitting region. In contrast to the previous fit, 
we observe here a behavior of the resulting cross section 
"■(piab) much more in line with our earlier analysis 5 . As 
a matter of fact we do not see any overprcdiction of the 
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FIG. 9: The result of the model employing fit strategy (II) for 
the real and imaginary part of the s-wave scattering amplitude 
of the r/N (top) and n~p — > rjn (bottom) channels. The error 
bands express the uncertainty due to variation of the model 
parameters as described in the text. 



total cross section at the position of the iV*(1535) peak 
which has been pointed out before, relying on the analy- 
sis of [13. 

As a further prediction we extract the scattering 
lengths of the r]N channels, which read 



1> V N- 



¥r)N = (+0.378™ + i0.20li°;°*?) f m , 



( 



2+0.092 
.101 

-o.208iS:21? 



-0.036^ 

■n 1OO+0.025 
iu - lo °-0.029 



) fm. (15) 



5 That means that no forced shift of the TV* (1535) pole to higher 
energies emerges from the fits. 



The observation to be made is that both are consistent 
with the extraction from other (more phenomcnological) 
approaches [39| . The s-wave amplitude in both channels 
can be found in Fig. [SJ For energies lower than the KA 
threshold we observe similar behavior as for the ampli- 
tudes extracted in fit strategy I, see Fig. [SJ 



F. Summary: Hadronic sector 

As an intermediate summary we conclude that the 
present framework serves as an appropriate theoretical 
tool for the analysis of hadronic scattering in a fairly 
wide energy range. Both isospin 1/2 s-wave resonances, 
the N*(1535) and A r *(1650) are generated dynamically 
from the NLO chiral potential in a quite natural way, 
i.e. resumming the Feynman diagrams in the BSE with 
the full off-shell dependence. At the same time it de- 
scribes the pion induced rj production of the nucleon in 
full agreement with the experimental data by Prakhov et 
al.. 
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FIG. 10: Types of diagrams of the turtle approximation. Shaded circles, squares and black dots represent the dressed meson- 
baryon vertex T, scattering amplitude T and the photon vertex as described in the main text. 



III. PION PHOTOPRODUCTION OFF THE 
PROTON 

Pion photo- and clcctroproduction off protons has been 
one of the premier objects of study in hadron physics for 
decades. One of the major issues was |U |42| and still 
is the interplay between the hadronic scattering and the 
photoproduction of mesons off the baryons. It is unques- 
tionable that the meson-baryon interaction plays a cru- 
cial role in the photoproduction processes via the rescat- 
tcring processes. At the production threshold, pion pho- 
toproduction can be successfully analysed within strictly 
perturbative ChPT as has been done to one loop about 
twenty years ago [H, |44|. Going to higher energies, 
one is again confronted with the problems already ap- 
pearing in the hadronic sector as discussed in the pre- 
vious section, namely, resonance phenomena. Thus a 
non-perturbative framework is required to implement the 
rescattering mechanism adequately. In the early years 
the unitarized hadronic amplitude was simply used as 
the final state interaction (FSI) multiplied on top of the 
7P — > ttN contact interaction, which in general violates 
the Ward-Takahashi identities and thus gauge invariance. 
Recently a framework for pion photoproduction based 
on the Julich model was constructed in [45[. There the 
hadronic part of the amplitude is also used as the FSI 
coupled to a special form of contact term and fulfilling 
the gauge Ward-Takahashi identities by construction. 

Our approach follows a different direction, where gauge 
invariance is not enforced via ad hoc conditions on ver- 
tex functions and propagators, but follows most straight- 
forwardly from the selected infinite subset of Feynman 
graphs which are summed up. The basic ideas can be 
traced back from Refs. [4ftf49| . Essentially, one adopts 
a generalization of the construction of a gauge-invariant 



amplitude (as spelled out e.g. in sec. (7.4) of |22|) to the 
present non-perturbative setting. It was first applied to 
the analysis of kaon photoproduction, relying on the lead- 
in g s-w ave terms from the three-flavor chiral Lagrangian, 
in [ljj] . There, the selected subset of Feynman graphs was 
referred to as "turtle approximation" . In principle this is 
the most natural way of constructing a gauge invariant 
photoproduction amplitude as the photon is coupled to 
any point of the p — > 4>B amplitude, ensuring current 
conservation. Nevertheless it requires as an input the 
underlying hadronic amplitude with the full off-shell de- 
pendence. Such an amplitude is provided in the previous 
section. It fulfills the two-body unitarity requirement ex- 
actly and the parameters are fixed such as to reproduce 
the s-wave of irN as well as irN — > rjN scattering. With- 
out any further fitting we wish to investigate what we can 
learn about the multipole amplitudes by just plugging in 
our fixed hadronic amplitude as an effective vertex func- 
tion. To put it in the words of Bcrcnds et al. [42j "we 
wish to see how the resonances from the hadronic spec- 
trum manifest themselves in the pion photoproduction 
observables" . 



A. Formalism 

Closely following the formalism explained in the 
gauge invariant photoproduction amplitude A^ M (q', k;p), 
is a sum of 8 different types of Feynman diagrams, see 
Fig. 1101 Here, q' is the four-momentum of the produced 
meson and p is again the overall four-momentum. The 
four-momentum of the incoming photon is denoted by k. 
The scattering amplitude T has been calculated in the 
previous section, consequently there are only two build- 
ings blocks left to be clarified, i.e. T denoting the dressed 
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FIG. 11: Symbolical representation of the dressed meson- 
baryon amplitude Y (circle), where the shaded square rep- 
resents the meson-baryon scattering amplitude T. 

meson-baryon vertex and the photon vertices W 7 0_>0, 

W^B-tB, W^B-^Bip and W^b^b^- 

The exact two-body unitarity is a crucial property of 
the hadronic amplitude. For this to be preserved in 
the photoproduction amplitude as well, the axial meson- 
baryon coupling has to be dressed properly. The tree 
level axial meson-baryon potential stems from the lead- 
ing order chiral Lagrangian Eq. ^) and reads 

V ax (q') = A x / 7 5 , (16) 

where we have separated off the channel space structure, 
which is specified in App. |XJ Dressing of this amplitude 
in the 'turtle approximation' [loj is presented in a rather 
intuitive pictorial way in Fig. QT] and reads 

T(q',p) =V ax {q') (17) 
+ / -^T(q>J-,p)iS(p- t)A(l)V ax (l), 

where p denotes the total four-momentum of this process 
which can take values of the proton momentum or the 
overall four-momentum of the photoproduction process. 
The scattering amplitude T consists of 20 different Dirac 
structures as presented in App. [B] and gives rise to 6 
different structures of the amplitude T, i.e. 

T( q ',p) = r x (p) • |5 7 5 + r a (p) • 7 5 + r 3 (p) • 4'i>i b 
+ r 4 (p)-d'i 5 + r 5 (p)- 1 5 (q' -p) 

+ T 6 (p)-ih 5 {q'-p). (18) 

The coefficients Tj(p) are elementary functions of p 2 , 
masses, scalar loop integrals, Im and Imb(P 2 ), collected 
in App \C\ coefficients Xj of the scattering amplitude as 
well as of the axial coupling constants D and F from 
cqn. ©. 

The numerical values of the latter constants should be 
taken with a grain of salt. Working at tree level for the 
axial vertex, the sum of D and F is set equal to the value 
of nuclcon axial vector charge g A , which has been in the 
focus of many experiments for decades and was measured 
recently to high precision in the neutron /3-decay using ul- 
tracold neutrons [H to be g A = 1.27590+ o ;^5- The ra- 
tio F/D is predicted by the SU(6) / non-relativistic quark 
model (NRQM) to be 2/3 which is actually quite close 
to the value of this ratio extracted from experiments, 
namely 0.58 ± 0.05, see [Fl| for a more detailed discus- 
sion. Within our approximation the axial coupling enters 



the photoproduction amplitudes dressed by the meson- 
baryon loops, for instance the kaon-loops. Those effects 
are known to be quite sizable and thus one has to choose 
at which level, i.e. on tree level or that of the dressed ver- 
tex r, one wishes to obtain an agreement with the physi- 
cal (measured) quantities. It turns out that although the 
value of r depends strongly on the choice of axial cou- 
plings, the photoproduction amplitudes calculated with 
both sets of axial couplings agree with each other within 
the uncertainty band. For this reason we stick to the 
commonly used values of D = 0.8 and F = 0.5. 

It remains now to specify how to couple the photon 
to the hadronic skeleton described above. For consis- 
tency reasons we shall consider the photon induced con- 
tact terms up to the second chiral order utilizing the first 
and second order chiral Lagrangian, Eq. ^ and Eq. ([7J . 
Previously we have set all the external currents to zero 
except the scalar one, whereas here we consider a vector 
current = — eQA^ with the electromagnetic vector po- 
tential and charge matrix Q = diag(2/3, —1/3,-1/3). 
A vector current modifies the covariant derivative as well 
as the chiral vielbein 

[D„,B] =d„B + -[([«t, d„u] - i{u\^u + uv^)),B], 
M/i =m t (9 AI C7 - i[Vfj,, U])v}. 

(2) 

A non- vanishing vector potential also features in C^g via 
the field-strength tensor 

/£" = u{d fl v l/ - d u v^ + u^(dV - d v v»)u. 

It appears in Eq. (JT)) accompanied by the two LECs &12 
and 613, which can not be determined from the scattering 
process. As mentioned before we do not wish to perform 
a fit for the photoproduction observables, therefore we 
stick to the values determined in Ref. [HH for these two 
new LECs. 

Furthermore, the pure mesonic chiral Lagrangian of 
second chiral order, which reads 

gives rise to the photon vertex W 7 ^_>^,. This and the 
remaining vertices are collected in App. lEI . 

Having specified all building blocks of the graphs col- 
lected in Fig. [TU] we now calculate the photoproduc- 
tion amplitude M.^ = Ylf=A^i' wnere the amplitudes 
Si correspond to a respective class of graphs defined in 
Fig. |T0j We wish to emphasize that there are 5 unitar- 
ity classes which by themselves obey two-body unitarity: 
{Sa] Sb+Se] Sc+Sd, Sf; Sq+Sh}- Gauge invariance is 
fulfilled for the amplitudes proportional to 612 and 613 au- 
tomatically. On the other hand for the remaining terms 
it is only fulfilled if all graphs presented in Fig. [TU] are 
taken into account, i.e. the photon is coupled to every 
possible part of the hadronic skeleton. 
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FIG. 12: Prediction for the multipole Eo+ for pion photoproduction corresponding to the hadronic solution I. For comparison, 
best fits of the MAID (circles) [531 an d SAID (squares) [5J| models are represented by blue and black symbols, respectively. 
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FIG. 13: Eq+ for the eta photoproduction as predicted based on hadronic solution I. For comparison, we also present the results 
of the ETAMAID (circles) [k| and Bonn-Gatchina (squares) analyses. 



B. Results 

In this section we present the prediction of the model 
for both sets of scattering amplitudes fixed in the 
hadronic sector in the last chapter. There in both strate- 
gies we have been concentrating on the description of the 
s-wave, thus we shall stick to the prediction of the quan- 
tities connected to this particular partial wave. Such 
a quantity is the electric multipole E$ +1 which can be 
expressed in terms of the Chew, Goldberger, Low and 
Nambu (CGLN) amplitudes as presented in App. [F] 

After having fixed the hadronic part of this ampli- 
tude, the photoproduction amplitude contains only 4 new 
parameters, namely D, F, 612 and 613. The first two 
are fixed to the commonly used values of D = 0.8 and 



F = 0.5. The 'magnetic' LECs &12 and 613 shall be taken 
from Ref. (52|, where they have been adjusted to fit the 
experimental data on magnetic moments. Within the un- 
certainty due to the choice of too there, these LECs are 
given by (in units of GcV -1 ) 

613 = 0.32 ±0.06 and b 12 = 0.095 ± 0.015. 

In order to give an impression of the uncertainty we 
proceed as follows. First for a fixed energy W cms and for 
each hadronic solution which lies in the uncertainty band 
of the hadronic solution as presented in the last chapter 
we calculate the photoproduction multipoles as functions 
of the 4 new LECs. Then for fixed D and F and a large 
set of randomly distributed values for b\2 and 613 within 
the uncertainty range on these two LECs we obtain a pre- 
diction on the photoproduction multipoles at the chosen 
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FIG. 14: Prediction for the multipole Ep+ for the pion photoproduction corresponding to the hadronic solution II. For com- 
parison best fits of the MAID (circles) [k| and SAID (squares) [54| models are represented by blue and black points with 
errorbars, respectively. 
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FIG. 15: i?o+ for the eta photoproduction as predicted starting from the hadronic solution II. For comparison we also present 
the outcome of the ETAMAID (circles) [k| and Bonn-Gatchina (squares) analysis. 



energy. Repeating this procedure for different energy val- 
ues we obtain a family of curves Eo+(W cms ). The hull of 
all those curves is assumed to reflect the uncertainty of 
the model properly. 

Solution I: As a central prediction of the present work 
we present the outcome of the multipole Eq + for pion 
photoproduction in the isospin 1/2 channel in Fig. 1121 
We restrict ourselves to energies below the second nu- 
clconic resonance as its position seems to be shifted as 
discussed in the last chapter. Without any fitting we 
observe an astonishing agreement of our prediction with 
the outcome of the fit from MAID2007 (updated uni tary 
isobar model) [53| and the one by the SAID group 54 1. 
Of course it is clear that, if none but irN channels are 
opened, Watson's theorem guarantees that the phase of 



E 0+ comes right, once the phase of ttN scattering has 
been fixed to the physical value. This theorem, however, 
does not fix the magnitude of the real and imaginary part 
of the photoproduction amplitude, neither it is clear how 
to apply it above the r/N threshold. 

The value of E$ + at the threshold has been debated 
for a long time, see [13 for a nice review on that topic. 
We obtain the following value 

K+{Sn) = (+10.4 ± 1.3) x 10- 3 /AU+, 

which has to be compared with E 0+ (Sn) = (+12.5 ± 
0.3) x 10~ 3 /A/ 7r + from experimental results |56l - l58| for a 
respective isospin combination. Seemingly our prediction 
is slightly below the the experimental result. Throughout 
this work we have not discussed the isospin 3/2 channel 
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and refrain from giving a numerical value of Eq + . For 
completeness, let us note that in all solutions we observe 
that the absolute value of the E 0+ (£31 ) is underestimated 
significantly. We trace that discrepancy to the missing 
Born graphs which are known to be important in this 
channel. We will come back to this issue in a forthcoming 
work [H. 

We can go further and make a prediction on the mul- 
tipolc amplitude for eta photoproduction. At the r]N 
threshold we extract the following value 



E v - 



(3.9±2.5) + i(10.7±2.7)) x lCn 3 / A 4+- 



The energy-dependence of Eo + (W cms ) is presented in 
Fig . IT31 where it is compared with fits by the ETAMAID 
|59( and Bonn-Gatchina JHO] groups. Seemingly there 
is a large qualitative agreement between our prediction 
and the phenomenological analysis by the ETAMAID 
and Bonn-Gatchina group. On a quantitative level we 
observe that the real part of the Eq + is suppressed com- 
pared to the outcome of the phenomenological analysis. 
We wish to remind the reader that the 'magnetic' LECs 
are taken from a tree level calculation only. In some addi- 
tional fits we have observed that the results of ETAMAID 
and Bonn-Gatchina group can be reproduced nicely in 
our approach using these LECs as free parameters. This, 
however, is not the original purpose of this work, namely 
the parameter free prediction of the photoproduction af- 
ter fixing the hadronic scattering. 

Solution II: Starting from the second hadronic solu- 
tion we obtain a prediction for pion photoproduction in 
the S11 channel as presented in Fig. [TU Although all 
parameters of the model are fixed in the hadronic solu- 
tion or taken from the literature as described above one 
observes a nice qualitative agreement of our prediction 
with the outcome of the SAID and MAID2007 analysis 
in a very large energy region. At the threshold we extract 
the following values for the lowest multipolc 

Eo + (S n ) = (+13.1 ±0.7) x lQ-yM^, 

which agrees nicely with the extraction from the experi- 
mental results in |56h58| . For higher energies, i.e. around 
1200 MeV and 1550 MeV, we observe a discrepancy of 
Eq+ compared to the fits by SAID and MAID groups. 
Additionally the uncertainty band appears quite under- 
estimated in this solution. We wish to remind the reader 
that the main difference of both solutions are the three 
regularization scales. In the first solution those are used 
as free fit parameters whereas in the second they are 
fixed. The particular choice of these values is motivated 
as described in the previous chapter, however, one should 
in principle investigate the influence of this choice on the 
hadronic solution as well as on the photoproduction am- 
plitudes. To do so one would have to refit the hadronic 
scattering for any other choice of the parameters in the 
solution II. Due to an enormous amount of computational 
time required for each fit, we refrain from including that 
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FIG. 16: A typical fit of our model to the SAID [54| (blue) 
and MAID [53| (black) analysis for the pion photoproduction 
as described in the text. 



uncertainty. One should keep in mind that a more re- 
alistic uncertainty band might be larger than the one 
presented here. 

To be complete we wish to comment now on the higher 
energy region, i.e. above the KA threshold, where the 
outcome of our prediction starts to deviate from the re- 
sults of the SAID and MAID groups. In fact this ob- 
servation is identical to the one made in the analysis of 
the photoproduction amplitudes in the Julich model [4(| , 
where no good overall fit could be achieved for the Eq + 
in the low and resonance energy region simultaneously. 
Although no fit to the photoproduction data was done 
in in the present work, four new parameters are enter- 
ing the calculation. The axial coupling as well as the 
'magnetic' LECs 612 and 613 are taken from estimations 
which rely on a strict perturbative calculation. Our non- 
perturbative framework is on the other hand suited to 
extend the range of applicability of the effective field the- 
ory. Thus it is a priori not clear whether it is sufficient 
to use this new LECs in the whole energy range. To 
underline this we fit our model to the SAID pion photo- 
production data with axial coupling and 'magnetic' LECs 
treated as free parameters. The best fit is presented in 
Fig. [TH1 where we observe a nice agreement above the KA 
threshold with the phenomenological models from SAID 
and MAID. 

For the eta photoproduction the prediction of the sec- 
ond solution is presented in Fig. \T5\ At the r/N threshold 
we obtain the following value 

E% + = f(-1.2 ±2.2) -H(6.9±2.3)) x lO^/A/^, 

which undershoots the numerical value obtained in the 
previous solution for the real part slightly and agrees for 
the imaginary part within the uncertainty range. The 
functional behavior of Eq+ is suppressed compared to 
the previous solution and even more in comparison to 
the ETAMAID and Bonn-Gatchina fits. As already dis- 
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cussed in the previous solution we can trace this discrep- 
ancy to the 'magnetic' LECs, which are taken from the 
tree level calculation [52| • These LECs do not change the 
functional form of the photoproduction amplitude but 
seem to enhance or suppress the structures present in 
the photoproduction amplitude. Those structures on the 
other hand seem to reflect one-to-one the structures aris- 
ing from the dynamics of the hadronic scattering process. 
Thus the correct description of meson photoproduction 
is necessarily to be connected to a proper description of 
the underlying hadronic scattering reactions. 
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Appendix A: Channel space structures 



For the channel indices {&, j;i,a} corresponding to the process 4>iB a — > 4>jBb^ the relevant coupling matrices from 
the leading, Eq. ([2]) , and next-to leading, Eq. ([7]) , order chiral Lagrangian read 



Ab,j;i,a 



A 



b,j;i,a 
14 



^(^[[A^nA ]), 

2 



A 



b,j;i,a 
ill 



bi ({X b Hx 3 \ [X\ \ a }}) + (A fc t[A\ [A jt , A a ]])) + fo 2 ((A bt {A jt , [A*, A a ]}) + (A bt {A\ [A* A a ]}) 

+ & 3 ((A bt {A^,{V,A Q }}) + (A 6t {A',{A^,A a }}>) + 26 4 (A bt A a )(A^A 
fe 5 (A 6t [[A jt , X% A a ]) + &6<A b t{[A^, A*], A a }) 

+ 6 7 ((A bt A 3t )(A 4 A a ) - (A 6t A')(A a A jt ) 
fe 8 ((A b t[A^'t ,[A J , A a ]]) + (A bt [A\ [A>'t, A a ]])) + & 9 ((A 6t [A^ ,{A' ; , A a }]) + (A b t[A l , {X j \ X a } 



Ab,j;i,a 
M 



2F j F i 



2b 



+ 6io((A 6 t{AJ't, {A 1 , A a }}) + (A fc t{A i , {A* A a }})J + 2b n (X b ^X a )(X^X 
. Q ((X^X a )([X^X l ]M)) + b D (<A b t{{A-''t, {M, A 1 }}, A a }> + (A bt {{A\ {M, A*}}, A Q })) 

+ b F ( (A b t [{A* {M,X 1 }}, X a ]) + (A b t [{A*, {M, A jt }}, X a ]) 



where A " denote the 3x3 channel matrices (e.g. <f> = <fr l X l for the physical meson fields) and the Fj are the meson 
decay constants in the respective channel. Moreover, M. is obtained from the quark mass matrix M. via the Gell-Mann 
Oakes Renner relations, and given in terms of the meson masses as follows, M = |diag(M|- + — M^-o + M^ ,Mf c0 — 
M 2 K+ + Ml ,M 2 K+ + M 2 K0 - M* ) . 

For the channel indices {&, j; a} corresponding to the process B a — > <pjB{, the channel-space matrix is given by 



.4 



b,j;a 



D 



(A bt {A 3t ,A a }) 



F 



V2Fj 



(A b [A jt , A a ]). 



Appendix B: Solution of the BSE with full off-shell dependence 

Here we present a technique for a solution of the Bethe Salpeter equation ([1]) with the full off-shell dependence. 
This method does not rely on any approximation of the BSE which are used very often in the literature, i.e. on-shell 
approximation or a three-dimensional reduction of the BSE to the Lippmann-Schwinger equation. It also is applicable 
for any kernel with only one restriction: The interaction kernel must consist of local terms only. Thus the solution of 
the BSE corresponds to an infinite chain of the Feynman bubble diagrams as presented in Fig. [TJ To keep this section 
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short wc will restrict the form of the kernel to the one used in the main text of this work. Up to the next-to-leading 
chiral order the meson-baryon local potential is given by the following form 

6 

V(q 2 ,qi;p) = J2 A i' D i(^,Qi;p), (Bl) 

i=l 

with 

v(q2,qi;p) = 92, (91-92), qkqA, t, {qk{q\-p), qk^i-p)), 

A = (A WT , A WT , {A 14 , + 2A 57 ), A 57 , A M , A SU ,A S11 ) , 

where and p denote the four momentum of the in-/outgoing meson and the overall four momentum of the 
scattering process. The capital letters A denote the channel space structures specified in the App. [A] where also the 
1 corresponds to a unity matrix. 

As already discussed the solution of the BSE corresponds to an infinite chain of Feynman diagrams which reads 

r fidl irk — f -|_ m 

T{q 2 ,q 1]P ) = V{q 2 ,q 1 ;p) + iJ — V { i2 J-p) {{ ^ i)2 _ m2){i2 ^ M2) V{l 4v p) + ... , 

where we again have suppressed the channel indices keeping in mind that T and V are matrices in channel space. From 
the last equation one easily sees that the iterative use of the interaction potential introduces new Dirac-momentum 
structures additionally to those of T>. The number of those is limited ensuring that one can separate the channel 
space structures from the Dirac-momentum structures of the scattering matrix as follows 

20 

T{q2,qv,p) = ^2 T?; ( s ) N i(?2,?i;p) , 
i=i 

with 

N(<?2,?i;p) p<fi, (ftftfii qfiqfi, fq/ifa-p), qk{q%-p), qh{q\-p), (fiq/i, {qi • p){qz • p), 

f{q\ - p){q2 - p), {qx-p),f{qi-p), (92 • qi),f{q2 ■ qi), qkp, qh, f{q2-p), fa-p), 1, f)- 

Please note that different to the decomposition of the potential V the coefficients Tj(s) still depend on the center-of- 
mass energy squared, i.e. s. The Dirac-momentum space spanned by the vectors I? is a subspace of those spanned 
by the vectors N. Consequently, the interaction potential V can also be expressed in terms of these vectors. More 
importantly 

f d d l a(p- I + m)b 

V« e *{q 2 ,l;p),be : J {27r)d {{p _ l)2 _ m2){p _ AP) e H^p). 

This ensures that Eq. ([T]) can be rewritten in a linear equation of the form 

X i j {s)T j {s) = V i for i,j = 1, .., 20 and X e Mat^ qi ' p)x ^' ll ' p l (B2) 

We refrain from presenting the matrix X as well as the solution of Eq. (|B2p explicitly for reasons of brevity but it 
should be mentioned that every element is a function of the baryon/meson masses, the cm. energy squared as well 
as of the scalar loop integrals Im and Imb, which are collected in the App. [Cj Moreover every element of X is a 
matrix in the channel space which depends on the LECs as presented in App. [5J The latter implies a non-commuting 
character of matrix elements X' , which one should keep in mind while solving Eq. (|B2[) . 



Appendix C: Loop integrals 

Here we collect all loop integrals required for the calculation of the scattering as well as the photoproduction 
amplitudes. Note that for reasons given in the main part all purely baryonic integrals are set to zero from the 
beginning. Utilizing dimensional regularization in the A/5 scheme the one-meson integral is given by 



Im{M) := 



d d l 



d=4 



1 



(27r) d I 2 - M 2 + k 



16tt 2 



2M 2 log(— ) 



MS 
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where fi is the regularization scale and M denotes the meson mass. We use in the following the common abbreviation 
A(a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 2bc, such that the meson-baryon (of masses M and m, respectively) integral 
reads 



Imb{s,M, m) 



d d l 1 



(2n) d I 2 - M 2 +ie(l- p) 2 - m 2 + ie 



MS 
d=4 1 



16tt 2 



i_i_oi t m \ i M 2 -m 2 + s ^ ,Ms y/\(s, M, m) ( y/X(s,M,m) ^ 

l + 21og(— ) + log(— ) -2-Z arctanhl v 

ji s m s V (m + My — s / 



The photoproduction amplitude involves further loop integrals. The triangle graph of class "D" in the Fig. [TUl gives 
rise to a meson-meson-baryon as well as via the Passarino-Veltman reduction to a meson-meson loop integral at 
s = k 2 , which read 



I MM {k 2 ,M) := 



d d l 



{2ir) d I 2 - M 2 +ie{l- k) 2 - M 2 + 



MS 



16tt 2 

lMMB(s,k 2 ,M,m) := 



- 1 r _ 2 + 2 log ( j£) - 2^^arctanh^V^™ 



i i 



(27r) d I 2 — M 2 + ie (I - k) 2 — M 2 +ie(l — k — pi) 2 - m 2 + ie ' 



MS 



where p\ = p — q' is the four-momentum of the incoming proton. We wish to emphasize that the photon coupled to 
a meson propagator does not induce a transition of this meson. Differently, coupled to a baryon propagator it can 
induce the £ •<-> A transition. Thus a meson-baryon-baryon loop integral required for the calculation of Feynman 
diagrams of class "E" reads in general 



lMBB(s,k ,M,mi,m 2 ) := 

MS 



d d l 



(27r) d I 2 - M 2 + ie (I - p{) 2 - m\ + ie (I - k - px) 2 - m| + ie 



Both last integrals can not be written in terms of elementary functions. We solve both integrals utilizing Cutkosky 
rules to calculate the imaginary part of the meson-baryon-baryon integral. A non-subtracted dispersion relation then 
gives the real part of the loop integral as follows, 

j fj , , 2 A , \\ 1 l ( m l - m 2 + k ' 2 - 2k 0l0 - 2g C ms&cms\ 

lm{I MBB (s,k ,M,m 1 ,m 2 )) = —— T log I — TTTi m i ' 

327rfc cms - v /s Vmf - + k A - 2k q + 2q cms k cms / 



1 f°° 

Re(lMBB(s,k 2 ,M, mi, 7112)) = — / ds 

77 J{m 2 +M) 2 



i lm(I MBB (s', k 2 ,M,m 1 ,m 2 )) 



^/A(s,pf,fc 2 )/(2yi), q cms = v / X(s,mlM 2 )/(2 y G) and fc = ^fc c 2 ms + k 2 , q = y/q 2 cms + M 2 . The same 
holds for the Immb, where in the last formulas one has to replace: mi — > M and mi — \ m. 



Appendix D: Partial wave analysis of nN — > rjN scattering 

The pion-induccd eta production off the neutron is dominated by the contribution of the nearby nucleon resonances, 
i.e. the <5ii(1535) and Di3(1520). From the previous study of elastic irN scattering we already know that the first 
one is described perfectly as a dynamically generated resonance within our approach, whereas the d-wave resonance 
is not. Thus we wish to clarify, whether an ansatz for the scattering amplitudes, which contains s- and p- waves only, 
is capable to generate a cos 2 (6>)-likc behavior of the differential cross sections da/dfl(s,cos(9)) (8 here denotes the 
scattering angle in the cm. frame). As a matter of fact this generation of a cos 2 (6>)-like structure through the iteration 
of p- waves does not seem to be appreciated in several experimental works, see e.g. [l4[. There, the presence of a 
cos 2 (#)-behavior in the shape of differential cross section is assumed to be a direct indication for a d-wavc dominance. 

Let us start from the most general form of the T-matrix, which is invariant under Lorentz as well as parity 
transformations. For the scattering of a meson-baryon system from initial state (i) with the meson momentum (<?), 
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and baryon momentum and spin (p, s) to the final state (/) with meson momentum (q'), and baryon momentum and 
spin (p' , s') it reads with the usual conventions used by Holder [6l[ 

M fi = -^u f {p',s'){A fi (s,t) + l($ + d')B fi (s,t)}u i (p,s), 

where s = P 2 := (p + q) 2 = (p' + q') 2 and t = (q — q') 2 = (p — p') 2 arc the Mandclstam variables. The amplitudes A 
and B can be recombined to the scattering amplitude on the mass shell Ton{$' , 4\ P) as follows 

T N($',<};P)=T£ )N (s,z)+f>TZ )N (s,z) = A(s,t) + ^ + $ , )B(s,t). 

Here, z = cos((9) is the standard representation of the scattering angle. In fact, z is related to the Mandelstam t via 

t = M 2 + M 2 - 2^q 2 + M 2 ^q 2 +M 2 + 2q t q f z , 

where q^if is the modulus of the center of mass momentum of the in- and outgoing system respectively. Suppressing 
the kinematic variables for the moment the differential cross section in the center of mass system reads 



For Dirac spinors normalized such that Uf(p)ui(p) = 2mdfi, and suppressing for brevity the channel indices, the spin 
sum can be calculated in terms of Tq N and Tq N as follows 

\ IW. s'){T% y + fT^Uiip, s) 2 = coolT^I 2 + 2c ai Re{T 1 o * N T ON ) + c n \T^ 



2 

s,s' 



n\ 2 ' 



where 



coo = y s (S s + m2 ~ M i)( s + m } ~ M f) + 4s(mimf - zg/gO) > 
coi = rrii(s + m 2 — M 2 ) + m/(s + m 2 — M 2 ), 

en = - ((s + m 2 - M 2 )(s + m)~ M 2 f ) + 4s(m i m / + zq f qijj . 

The above formulae specify all required kinematics and spin structure. The dynamical input is incorporated within 
the scattering amplitudes, A and B. In the main body of this work these are taken to be solutions of the BSE. In view 
of the above question, we wish to make an ansatz for the scattering amplitudes. First of all, the standard amplitudes 
A and B can be expanded in Legendre polynomials Pi{z) as follows 



A(s,t) 




+ m t 


47T 


\J ^cms;j 


+ mi 


B(s,t) 


l 




4-7T 


\J E cm s;i 


+ Vtli 



V "cms m f 

fl{S,t)—= :J2(M) 





h mf yJWcrns - rrn 


\J E cms -f 


+ m f V Pcms;i - rrii 


1 


1 

1 



7= — + f 2 (a,t) 



1 



y/E cma -J + TUf \/E cmS ;i mi \/E cras -f irif 



where Wcms = yfs and -B cms: i// = yqf/f + m 1/f After a variable transformation the amplitudes f± t 2 are related to 
the commonly used partial wave amplitudes fi±(s) as follows (6ll. |62|. 

oo 

Ms, Z) = £(/(,_!)+(«) - f (l+1) -( S ))Pl(z), 

1 = 1 

OO 

f 2 (s,z) = J2(fl-(s)-fi+(s))P!(z)- 
i=i 

For the purpose of this section we do not consider additional constraints for the partial wave amplitudes, e.g. due 
to analyticity or unitarity. Thus both real and imaginary part of those are used as free parameters, which will be 
adjusted to reproduce the data on differential cross sections for the process ir~N r/N, measured by Prakhov et al., 
see Ref. [Tij . For the truncation of the partial wave expansion we assume three different scenarios: 
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FIG. 17: Best fits of three assumed scenarios to the data of the differential cross sections from Ref. [3] for different pion momenta 
piab- The dashed (blue), red (dotted) and green (full) line correspond to the first, second and third scenario, respectively, as 
described in the text. 



1. The scattering amplitude contains only the s-wave dynamically, meant in the sense of the above discussion. Free 
parameters are {Re(/o+), Im(/ +)}. 

2. Both, s- and p- waves are included dynamically. This is the case for the solution of the BSE with contact terms 
from NLO chiral Lagrangian as performed in our approach. Without restricting the parity of the p- waves we end 
up with the following free parameters for this scenario: {Re(/o+), Im(/o+), Re(/i _), Im(/i_), Re(/i_|_), Im(/i + )}. 

3. The scattering amplitude is determined by s- and d-wave (for instance -D13), whereas the p-wave is negli- 
gible. This is the case for the process irN — > rjN from the phenomenological point of view. For phe- 
nomenological reasons we assume only the D13 wave to be non- negligible. Thus the free parameters are 
{Re(/o + ),Im(/o + ),Re(/ 2 _),Im(/ 2 _)}. 

For each incident 7r _ momenta separately and for each scenario we obtain best fits as presented in Fig. [T7] As 
expected the first scenario is only capable to fit the data at lowest beam momenta. The s-wave is dominant at low 
energies, however at higher energies it lacks the angular dependence and thus fails to describe the data properly. 
Going to higher beam momenta both, the second and third scenario describe the data equally well. It turns out that 
the presence of p-waves of both parities is required to reproduce the z 2 -behavior. Thus albeit our approach based 
on the unitarization of the NLO chiral potential does not produce d-waves in the sense of the above discussion, it is 
capable to reproduce the data on differential cross section well enough. 
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Appendix E: One-photon vertices 

For all vertices the in(out)going meson and baryon are denoted by the channel index i(J) and a(b) respectively, 
whereas the charge of the corresponding particle is denoted by Q. All required photon induced vertices W from the 
leading order chiral Lagrangian Eq. (O read 

where qi denotes the four-momentum of the produced (outgoing) meson. The vertices from the second order chiral 
Lagrangian (J7J posses more involved channel structures, which for instance can be traced back elegantly to the already 
defined channel matrices inlAl The interaction vertex for the process j(k)4>i(qi)B a (jpi) — > 4>j{q_2)Bb{p2) reads 

+ tA 57 - (Qi[y.^] + Oi[y.0]) 

+ i;4«ii • (Qj^(qi,pi + qi) + Q,qfi{q^ + + Q^fe, 92 + Pa) + Qi$(q£ +Pa) + #«i + #«a) 
- ^(6i2(A 6 t[Q4A\A^]-Q i [A^,A i ],Al) + 6i 3 (A 6 t{Q4A\A^]-Q j [A^^A i ],A a })) ■ (W-V*). 

The latter expression originates from the electromagnetic field-strength tensor fl?" '. Furthermore the same term gives 
rise to an additional coupling of the photon to a baryon, which docs not vanish for electrically neutral baryons. It 
also induces a baryon transition £ A, the corresponding vertex reads 

Kb^b = 2*e(b 12 (^[Q,\ a }) +b 13 (X b HQ,X a })) ■ (jfry"-7 M *), 
where Q = diag(2/3, —1/3, —1/3) is the charge matrix and e is the charge of an electron. 



Appendix F: Multipoles 

In this section we wish to specify the major technical steps on the way from the photoproduction amplitude 
as calculated utilizing usual Feynman rules to the multipole amplitudes as well as to the cross sections. In large 
parts of this section we use the conventions of |42| and start from the most general Lorentz covariant transition 
matrix element for the process of meson (<pf) production of the baryon (Bi) via an incoming photon (7(fc)), i.e. 
7(fc)i3 ? (p— k) — > Bf(p— q)<pf(q). It reads 

8 

Tf^ie^UfiY^B^m, (Fl) 
fc=i 

where e M is the photon polarization vector. The initial and final Dirac spinors Ui and Uf are normalized like 
uu = 2m, with m the mass of the corresponding baryon. The coefficients Bi are functions of the coefficients 
of the hadronic scattering amplitude {T^;i = 1,...,20} as defined in the App. [Bj loop integrals from App. [C] 
as well as vertices from App. [K\ Since both baryons arc on-shcll there are only 8 different structures, i.e. 
K G {757^, 27^,275^,275^, 757^,75^,75^, 75^} withP=i(2p-g-fe). 

Fixing the axis of quantization to the z-axis, one is able to reduce the Dirac spinors to the two-component spinors 
X as follows 

8 

T fi = 8ir^ x) Yl F*8kXi ■ (F2) 
fc=i 

This gives rise to the so-called Chew, Goldbcrgcr, Low and Nambu (CGLN) [62j amplitudes Fi, which are defined 
in the basis given by Q k <G ■ •q)(& '[kx e]),i(<r -k)(q- e),i(<7 -q)(q- ~e),i((T -k)(k- e),i(a -q)(k- ~t 
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),i(cr -q){k- e),i(cr -g)eo,i(cr -fc)eo}. Here, an arrow denotes a three-dimensional vector and a hat a normalized 
three-vector. Due to current conservation, two of the eight CGLN amplitudes can be eliminated via 

+ q)Tz + J5 - jr^n = and (k ■ q)F 4 + Fq - j-: T7 = 0, 
\k\ \k\ 

which serves as a good check of our calculation. To further extent two of the remaining six amplitudes are accompanied 
by scalar components of e only and thus have no influence on photoproduction amplitudes, i.e. process including real 
photons. Finally the lowest electric multipolc Eo + can be calculated as follows 

E 0+ = J i dz^Po?! - +\i p o- W2) , (F3) 

where Pi denote the Legendre polynomials. The latter as well as the CGLN amplitudes are functions of the cosine of 
the scattering angle in the cm. system, z. The unpolarized differential cross section for meson photoproduction is 
given by 

% = m(\^\ 2 + W + + \\?a\ 2 + Re(^) + Rc(J- 2 J" 3 *) + zRc(T 3 T: - 2T X T*) 

- z 2 {\\^? + \\H 2 + Rc(J"i^ 4 * + ^2^3)) - z 3 Re(J- 3 J|)) . (F4) 
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